
function [mout] = dsp(num)

% ****************************************************************
% num       : 输出的差集偶的个数
% mout      : 输出的差集偶
% ****************************************************************
% 差集偶，长度为13[[1, 3, 6, 7, 8, 9, 12], [3, 6, 7, 8, 9, 12]]
% 对应的特征序列偶[[1,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1],[-1,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1]]
mout = {[[1,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1];[-1,-1,1,-1,-1,1,1,1,1,-1,-1,1,-1]]};
if num > 1
    for ii=2:num
        mout{end+1} = [shift(mout{ii-1}(1,:), 1);shift(mout{ii-1}(2,:), 1)];
    end
end